Long-read transcriptome sequencing provides insight into lignan biosynthesis during fruit development in Schisandra chinensis.

BACKGROUND: Schisandra chinensis, an ancient member of the most basal angiosperm lineage which is known as the ANITA, is a fruit-bearing vine with the pharmacological effects of a multidrug system, such as antioxidant, anti-inflammatory, cardioprotective, neuroprotective, anti-osteoporosis effects. Its major bioactive compound is represented by lignans such as schisandrin. Molecular characterization of lignan biosynthesis in S. chinensis is of great importance for improving the production of this class of active compound. However, the biosynthetic mechanism of schisandrin remains largely unknown. RESULTS: To understand the potential key catalytic steps and their regulation of schisandrin biosynthesis, we generated genome-wide transcriptome data from three different tissues of S. chinensis cultivar Cheongsoon, including leaf, root, and fruit, via long- and short-read sequencing technologies. A total of 132,856 assembled transcripts were generated with an average length of 1.9 kb and high assembly completeness. Overall, our data presented effective, accurate gene annotation in the prediction of functional pathways. In particular, the annotation revealed the abundance of transcripts related to phenylpropanoid biosynthesis. Remarkably, transcriptome profiling during fruit development of S. chinensis cultivar Cheongsoon revealed that the phenylpropanoid biosynthetic pathway, specific to coniferyl alcohol biosynthesis, showed a tendency to be upregulated at the postfruit development stage. Further the analysis also revealed that the pathway forms a transcriptional network with fruit ripening-related genes, especially the ABA signaling-related pathway. Finally, candidate unigenes homologous to isoeugenol synthase 1 (IGS1) and dirigent-like protein (DIR), which are subsequently activated by phenylpropanoid biosynthesis and thus catalyze key upstream steps in schisandrin biosynthesis, were identified. Their expression was increased at the postfruit development stage, suggesting that they may be involved in the regulation of schisandrin biosynthesis in S. chinensis. CONCLUSIONS: Our results provide new insights into the production and accumulation of schisandrin in S. chinensis berries and will be utilized as a valuable transcriptomic resource for improving the schisandrin content.

Influence of Genotype on High Glucosinolate Synthesis Lines of Brassica rapa.

This study was conducted to investigate doubled haploid (DH) lines produced between high GSL (HGSL) Brassica rapa ssp. trilocularis (yellow sarson) and low GSL (LGSL) B. rapa ssp. chinensis (pak choi) parents. In total, 161 DH lines were generated. GSL content of HGSL DH lines ranged from 44.12 to 57.04 µmol·g-1·dry weight (dw), which is within the level of high GSL B. rapa ssp. trilocularis (47.46 to 59.56 µmol g-1 dw). We resequenced five of the HGSL DH lines and three of the LGSL DH lines. Recombination blocks were formed between the parental and DH lines with 108,328 single-nucleotide polymorphisms in all chromosomes. In the measured GSL, gluconapin occurred as the major substrate in HGSL DH lines. Among the HGSL DH lines, BrYSP_DH005 had glucoraphanin levels approximately 12-fold higher than those of the HGSL mother plant. The hydrolysis capacity of GSL was analyzed in HGSL DH lines with a Korean pak choi cultivar as a control. Bioactive compounds, such as 3-butenyl isothiocyanate, 4-pentenyl isothiocyanate, 2-phenethyl isothiocyanate, and sulforaphane, were present in the HGSL DH lines at 3-fold to 6.3-fold higher levels compared to the commercial cultivar. The selected HGSL DH lines, resequencing data, and SNP identification were utilized for genome-assisted selection to develop elite GSL-enriched cultivars and the industrial production of potential anti-cancerous metabolites such as gluconapin and glucoraphanin.

Development of 50 InDel-based barcode system for genetic identification of tartary buckwheat resources.

Tartary buckwheat (Fagopyrum tataricum Gartn.) is a highly functional crop that is poised to be the target of many future breeding efforts. The reliable ex situ conservation of various genetic resources is essential for the modern breeding of tartary buckwheat varieties. We developed PCR-based co-dominant insertion/deletion (InDel) markers to discriminate tartary buckwheat genetic resources. First, we obtained the whole genome from 26 accessions across a superscaffold-scale reference genome of 569.37 Mb for tartary buckwheat cv. "Daegwan 3-7." Next, 171,926 homogeneous and 53,755 heterogeneous InDels were detected by comparing 26 accessions with the "Daegwan 3-7" reference sequence. Of these, 100 candidate InDels ranging from 5-20 bp in length were chosen for validation, and 50 of them revealed polymorphisms between the 26 accessions and "Daegwan 3-7." The validated InDels were further tested through the assessment of their likelihood to give rise to a single or a few PCR products in 50 other accessions, covering most tartary buckwheat genome types. The major allele frequencies ranged from 0.5616 at the TB42 locus to 0.9863 at the TB48 locus, with the average PIC value of 0.1532 with a range of 0.0267-0.3712. To create a user-friendly system, the homology of the genotypes between and among the accessions were visualized in both one- (1D) and two-dimensional (2D) barcode types by comparing amplicon polymorphisms with the reference variety, "Daegwan 3-7." A phylogenetic tree and population structure of the 76 accessions according to amplicon polymorphisms for the 50 InDel markers corresponded to those using non-synonymous single nucleotide polymorphism variants, indicating that the barcode system based on the 50 InDels was a useful tool to improve the reliability of identification of tartary buckwheat accessions in the germplasm stocks.

Integrated genomic and transcriptomic analysis reveals unique mechanisms for high osmotolerance and halotolerance in Hyphopichia yeast.

The yeast species Hyphopichia is common in nature and strongly competitive under harsh environmental conditions. Here, we characterized Hyphopichia burtonii KJJ43 and H. pseudoburtonii KJS14, which exhibit strong halotolerance, using genomic and transcriptomic analyses. The genomes of H. burtonii and H. pseudoburtonii comprised eight chromosomes with 85.17% nucleotide identity and significant divergence in synteny. Notably, both Hyphopichia genomes possessed extended gene families of amino acid permeases and ATP-binding cassette (ABC) transporters, whose dynamic expression patterns during osmotic stress were revealed using transcriptome profiling. Intriguingly, we found unique features of the HOG pathway activated by Hog1p even under non-osmotic stress conditions and the upregulation of cytosolic Gpd1 protein during osmotic stress. Associated with hyperfilamentation growth under high osmotic conditions, a set of genes in the FLO family with induced expression in response to NaCl, KCl, and sorbitol supplementation were identified. Moreover, comparative transcriptome analysis reveals the NaCl-specific induction of genes involved in amino acid biosynthesis and metabolism, particularly BAT2. This suggests the potential association between oxoacid reaction involving branched-chain amino acids and osmotolerance. The combined omics analysis of two Hyphopichia species provides insights into the novel mechanisms involved in salt and osmo-stress tolerance exploited by diverse eukaryotic organisms.

Brassinosteroid-BZR1/2-WAT1 module determines the high level of auxin signalling in vascular cambium during wood formation.

The tight regulation of local auxin homeostasis and signalling maxima in xylem precursor cells specifies the organising activity of the vascular cambium and consequently promotes xylem differentiation and wood formation. However, the molecular mechanisms underlying the local auxin signalling maxima in the vascular cambium are largely unknown. Here, we reveal that brassinosteroid (BR)-activated WALLS ARE THIN1 (WAT1) facilitates wood formation by enhancing local auxin signalling in the vascular cambium in Solanum lycopersicum. Growth defects and low auxin signalling readouts in the BR-deficient tomato cultivar, Micro-Tom, were associated with a novel recessive allele, Slwat1-copi, created by the insertion of a retrotransposon in the last exon of the SlWAT1 locus. Molecular and genetic studies by generating the gain-of-function and loss-of-function tomato mutants revealed that SlWAT1 is a critical regulator for fine tuning local auxin homeostasis and signalling outputs in vascular cambium to facilitate secondary growth. Finally, we discovered that BR-regulated SlBZR1/2 directly activated downstream auxin responses by SlWAT1 upregulation in xylem precursor cells to facilitate xylem differentiation and subsequent wood formation. Our data suggest that the BR-SlBZR1/2-WAT1 signalling network contributes to the high level of auxin signalling in the vascular cambium for secondary growth.

Whole-genome, transcriptome, and methylome analyses provide insights into the evolution of platycoside biosynthesis in Platycodon grandiflorus, a medicinal plant.

Triterpenoid saponins (TSs) are common plant defense phytochemicals with potential pharmaceutical properties. Platycodon grandiflorus (Campanulaceae) has been traditionally used to treat bronchitis and asthma in East Asia. The oleanane-type TSs, platycosides, are a major component of the P. grandiflorus root extract. Recent studies show that platycosides exhibit anti-inflammatory, antiobesity, anticancer, antiviral, and antiallergy properties. However, the evolutionary history of platycoside biosynthesis genes remains unknown. In this study, we sequenced the genome of P. grandiflorus and investigated the genes involved in platycoside biosynthesis. The draft genome of P. grandiflorus is 680.1 Mb long and contains 40,017 protein-coding genes. Genomic analysis revealed that the CYP716 family genes play a major role in platycoside oxidation. The CYP716 gene family of P. grandiflorus was much larger than that of other Asterid species. Orthologous gene annotation also revealed the expansion of ß-amyrin synthases (bASs) in P. grandiflorus, which was confirmed by tissue-specific gene expression. In these expanded gene families, we identified key genes showing preferential expression in roots and association with platycoside biosynthesis. In addition, whole-genome bisulfite sequencing showed that CYP716 and bAS genes are hypomethylated in P. grandiflorus, suggesting that epigenetic modification of these two gene families affects platycoside biosynthesis. Thus whole-genome, transcriptome, and methylome data of P. grandiflorus provide novel insights into the regulation of platycoside biosynthesis by CYP716 and bAS gene families.

Development of 18 microsatellite markers for Atractylodes japonica.

PREMISE: Atractylodes japonica (Asteraceae) is endemic to East Asia, where its rhizomes are used in traditional medicine. To investigate the genetic diversity of this species, we developed polymorphic microsatellite markers. METHODS AND RESULTS: We obtained a total of 175,825 simple sequence repeat (SSR) loci using the Illumina HiSeq 2500 system. Eighteen polymorphic SSR primer pairs were selected to determine heterozygosity levels and allele numbers in 80 individuals from four A. japonica populations. The levels of observed and expected heterozygosity ranged from 0.000 to 1.000 and from 0.133 to 0.892, respectively. Cross-amplification in the related species A. macrocephala and A. lancea was successful in 15 and 14 of the 18 markers, respectively. CONCLUSIONS: These microsatellite markers will be useful for future studies involving A. japonica population genetics and breeding.

Efficient and specific generation of knockout mice using Campylobacter jejuni CRISPR/Cas9 system.

The Streptococcus pyogenes CRISPR/Cas9 (SpCas9) system is now widely utilized to generate genome engineered mice; however, some studies raised issues related to off-target mutations with this system. Herein, we utilized the Campylobacter jejuni Cas9 (CjCas9) system to generate knockout mice. We designed sgRNAs targeting mouse Tyr or Foxn1 and microinjected into zygotes along with CjCas9 mRNA. We obtained newborn mice from the microinjected embryos and confirmed that 50% (Tyr) and 38.5% (Foxn1) of the newborn mice have biallelic mutation on the intended target sequences, indicating efficient genome targeting by CjCas9. In addition, we analyzed off-target mutations in founder mutant mice by targeted deep sequencing and whole genome sequencing. Both analyses revealed no off-target mutations at potential off-target sites predicted in silico and no unexpected random mutations in analyzed founder animals. In conclusion, the CjCas9 system can be utilized to generate genome edited mice in a precise manner.

TaF: a web platform for taxonomic profile-based fungal gene prediction.

INTRODUCTION: The accurate prediction and annotation of gene structures from the genome sequence of an organism enable genome-wide functional analyses to obtain insight into the biological properties of an organism. OBJECTIVES: We recently developed a highly accurate filamentous fungal gene prediction pipeline and web platform called TaF. TaF is a homology-based gene predictor employing large-scale taxonomic profiling to search for close relatives in genome queries. METHODS: TaF pipeline consists of four processing steps; (1) taxonomic profiling to search for close relatives to query, (2) generation of hints for determining exon-intron boundaries from orthologous protein sequence data of the profiled species, (3) gene prediction by combination of ab inito and evidence-based prediction methods, and (4) homology search for gene models. RESULTS: TaF generates extrinsic evidence that suggests possible exon-intron boundaries based on orthologous protein sequence data, thus reducing false-positive predictions of gene structure based on distantly related orthologs data. In particular, the gene prediction method using taxonomic profiling shows very high accuracy, including high sensitivity and specificity for gene models, suggesting a new approach for homology-based gene prediction from newly sequenced or uncharacterized fungal genomes, with the potential to improve the quality of gene prediction. CONCLUSION: TaF will be a useful tool for fungal genome-wide analyses, including the identification of targeted genes associated with a trait, transcriptome profiling, comparative genomics, and evolutionary analysis.

Long-term health and germline transmission in transgenic cattle following transposon-mediated gene transfer.

BACKGROUND: Transposon-mediated, non-viral gene delivery is a powerful tool for generating stable cell lines and transgenic animals. However, as multi-copy insertion is the preferred integration pattern, there is the potential for uncontrolled changes in endogenous gene expression and detrimental effects in cells or animals. Our group has previously reported on the generation of several transgenic cattle by using microinjection of the Sleeping Beauty (SB) and PiggyBac (PB) transposons and seeks to explore the long-term effects of this technology on cattle. RESULTS: Transgenic cattle, one female (SNU-SB-1) and one male (SNU-PB-1), reached over 36 months of age with no significant health issues and normal blood parameters. The detection of transgene integration and fluorescent signal in oocytes and sperm suggested the capacity for germline transmission in both of the founder animals. After natural breeding, the founder transgenic cow delivered a male calf and secreted milk containing fluorescent transgenic proteins. The calf expressed green fluorescent protein in primary cells from ear skin, with no significant change in overall genomic stability and blood parameters. Three sites of transgene integration were identified by next-generation sequencing of the calf's genome. CONCLUSIONS: Overall, these data demonstrate that transposon-mediated transgenesis can be applied to cattle without being detrimental to their long-term genomic stability or general health. We further suggest that this technology may be usefully applied in other fields, such as the generation of transgenic animal models.

Development of Genome-Wide SSR Markers from Angelica gigas Nakai Using Next Generation Sequencing.

Angelica gigas Nakai is an important medicinal herb, widely utilized in Asian countries especially in Korea, Japan, and China. Although it is a vital medicinal herb, the lack of sequencing data and efficient molecular markers has limited the application of a genetic approach for horticultural improvements. Simple sequence repeats (SSRs) are universally accepted molecular markers for population structure study. In this study, we found over 130,000 SSRs, ranging from di- to deca-nucleotide motifs, using the genome sequence of Manchu variety (MV) of A. gigas, derived from next generation sequencing (NGS). From the putative SSR regions identified, a total of 16,496 primer sets were successfully designed. Among them, we selected 848 SSR markers that showed polymorphism from in silico analysis and contained tri- to hexa-nucleotide motifs. We tested 36 SSR primer sets for polymorphism in 16 A. gigas accessions. The average polymorphism information content (PIC) was 0.69; the average observed heterozygosity (HO) values, and the expected heterozygosity (HE) values were 0.53 and 0.73, respectively. These newly developed SSR markers would be useful tools for molecular genetics, genotype identification, genetic mapping, molecular breeding, and studying species relationships of the Angelica genus.

Isoform Sequencing Provides a More Comprehensive View of the Panax ginseng Transcriptome.

Korean ginseng (Panax ginseng C.A. Meyer) has been widely used for medicinal purposes and contains potent plant secondary metabolites, including ginsenosides. To obtain transcriptomic data that offers a more comprehensive view of functional genomics in P. ginseng, we generated genome-wide transcriptome data from four different P. ginseng tissues using PacBio isoform sequencing (Iso-Seq) technology. A total of 135,317 assembled transcripts were generated with an average length of 3.2 kb and high assembly completeness. Of those unigenes, 67.5% were predicted to be complete full-length (FL) open reading frames (ORFs) and exhibited a high gene annotation rate. Furthermore, we successfully identified unique full-length genes involved in triterpenoid saponin synthesis and plant hormonal signaling pathways, including auxin and cytokinin. Studies on the functional genomics of P. ginseng seedlings have confirmed the rapid upregulation of negative feed-back loops by auxin and cytokinin signaling cues. The conserved evolutionary mechanisms in the auxin and cytokinin canonical signaling pathways of P. ginseng are more complex than those in Arabidopsis thaliana. Our analysis also revealed a more detailed view of transcriptome-wide alternative isoforms for 88 genes. Finally, transposable elements (TEs) were also identified, suggesting transcriptional activity of TEs in P. ginseng. In conclusion, our results suggest that long-read, full-length or partial-unigene data with high-quality assemblies are invaluable resources as transcriptomic references in P. ginseng and can be used for comparative analyses in closely related medicinal plants.

accD nuclear transfer of Platycodon grandiflorum and the plastid of early Campanulaceae.

BACKGROUND: Campanulaceae species are known to have highly rearranged plastid genomes lacking the acetyl-CoA carboxylase (ACC) subunit D gene (accD), and instead have a nuclear (nr)-accD. Plastid genome information has been thought to depend on studies concerning Trachelium caeruleum and genome announcements for Adenophora remotiflora, Campanula takesimana, and Hanabusaya asiatica. RNA editing information for plastid genes is currently unavailable for Campanulaceae. To understand plastid genome evolution in Campanulaceae, we have sequenced and characterized the chloroplast (cp) genome and nr-accD of Platycodon grandiflorum, a basal member of Campanulaceae. RESULTS: We sequenced the 171,818 bp cp genome containing a 79,061 bp large single-copy (LSC) region, a 42,433 bp inverted repeat (IR) and a 7840 bp small single-copy (SSC) region, which represents the cp genome with the largest IR among species of Campanulaceae. The genome contains 110 genes and 18 introns, comprising 77 protein-coding genes, four RNA genes, 29 tRNA genes, 17 group II introns, and one group I intron. RNA editing of genes was detected in 18 sites of 14 protein-coding genes. Platycodon has an IR containing a 3' rps12 operon, which occurs in the middle of the LSC region in four other species of Campanulaceae (T. caeruleum, A. remotiflora, C. takesimana, and H. asiatica), but lacks accD, clpP, infA, and rpl23, as has been found in these four species. Platycodon nr-accD contains about 3.2 kb intron between nr-accD.e1 and nr-accD.e2 at the same insertion point as in other Campanulaceae. The phylogenies of the plastid genomes and accD show that Platycodon is basal in the Campanulaceae clade, indicating that IR disruption in Campanulaceae occurred after the loss of accD, clpP, infA, and rpl23 in the cp genome, which occurred during plastid evolution in Campanulaceae. CONCLUSIONS: The plastid genome of P. grandiflorum lacks the rearrangement of the IR found in T. caeruleum, A. remotiflora, C. takesimana, and H. asiatica. The absence of accD, clpP, infA, and rpl23 in the plastid genome is a synapomorphic characteristic of Campanulaceae. The chloroplast genome phylogeny supports the hypothesis that chloroplast genomic arrangement occurred after accD nuclear transfer and loss of the four genes in the plastid of early Campanulaceae as a lineage of asterids.

Long-read transcriptome data for improved gene prediction in Lentinula edodes.

Lentinula edodes is one of the most popular edible mushrooms in the world and contains useful medicinal components such as lentinan. The whole-genome sequence of L. edodes has been determined with the objective of discovering candidate genes associated with agronomic traits, but experimental verification of gene models with correction of gene prediction errors is lacking. To improve the accuracy of gene prediction, we produced 12.6 Gb of long-read transcriptome data of variable lengths using PacBio single-molecule real-time (SMRT) sequencing and generated 36,946 transcript clusters with an average length of 2.2 kb. Evidence-driven gene prediction on the basis of long- and short-read RNA sequencing data was performed; a total of 16,610 protein-coding genes were predicted with error correction. Of the predicted genes, 42.2% were verified to be covered by full-length transcript clusters. The raw reads have been deposited in the NCBI SRA database under accession number PRJNA396788.

Whole-genome de novo sequencing, combined with RNA-Seq analysis, reveals unique genome and physiological features of the amylolytic yeast Saccharomycopsis fibuligera and its interspecies hybrid.

BACKGROUND: Genomic studies on fungal species with hydrolytic activity have gained increased attention due to their great biotechnological potential for biomass-based biofuel production. The amylolytic yeast Saccharomycopsis fibuligera has served as a good source of enzymes and genes involved in saccharification. Despite its long history of use in food fermentation and bioethanol production, very little is known about the basic physiology and genomic features of S. fibuligera. RESULTS: We performed whole-genome (WG) de novo sequencing and complete assembly of S. fibuligera KJJ81 and KPH12, two isolates from wheat-based Nuruk in Korea. Intriguingly, the KJJ81 genome (~38 Mb) was revealed as a hybrid between the KPH12 genome (~18 Mb) and another unidentified genome sharing 88.1% nucleotide identity with the KPH12 genome. The seven chromosome pairs of KJJ81 subgenomes exhibit highly conserved synteny, indicating a very recent hybridization event. The phylogeny inferred from WG comparisons showed an early divergence of S. fibuligera before the separation of the CTG and Saccharomycetaceae clades in the subphylum Saccharomycotina. Reconstructed carbon and sulfur metabolic pathways, coupled with RNA-Seq analysis, suggested a marginal Crabtree effect under high glucose and activation of sulfur metabolism toward methionine biosynthesis under sulfur limitation in this yeast. Notably, the lack of sulfate assimilation genes in the S. fibuligera genome reflects a unique phenotype for Saccharomycopsis clades as natural sulfur auxotrophs. Extended gene families, including novel genes involved in saccharification and proteolysis, were identified. Moreover, comparative genome analysis of S. fibuligera ATCC 36309, an isolate from chalky rye bread in Germany, revealed that an interchromosomal translocation occurred in the KPH12 genome before the generation of the KJJ81 hybrid genome. CONCLUSIONS: The completely sequenced S. fibuligera genome with high-quality annotation and RNA-Seq analysis establishes an important foundation for functional inference of S. fibuligera in the degradation of fermentation mash. The gene inventory facilitates the discovery of new genes applicable to the production of novel valuable enzymes and chemicals. Moreover, as the first gapless genome assembly in the genus Saccharomycopsis including members with desirable traits for bioconversion, the unique genomic features of S. fibuligera and its hybrid will provide in-depth insights into fungal genome dynamics as evolutionary adaptation.

Whole genome de novo sequencing and genome annotation of the world popular cultivated edible mushroom, Lentinula edodes.

Lentinula edodes, the popular shiitake mushroom, is one of the most important cultivated edible mushrooms. It is used as a food and for medicinal purposes. Here, we present the 46.1 Mb draft genome of L. edodes, comprising 13,028 predicted gene models. The genome assembly consists of 31 scaffolds. Gene annotation provides key information about various signaling pathways and secondary metabolites. This genomic information should help establish the molecular genetic markers for MAS/MAB and increase our understanding of the genome structure and function.

Obesity Alters the Microbial Community Profile in Korean Adolescents.

Obesity is an increasing public health concern worldwide. According to the latest Organization for Economic Co-operation and Development (OECD) report (2014), the incidence of child obesity in Korea has exceeded the OECD average. To better understand and control this condition, the present study examined the composition of the gut microbial community in normal and obese adolescents. Fecal samples were collected from 67 obese (body mass index [BMI] ≥ 30 kg/m2, or ≥ 99th BMI percentile) and 67 normal (BMI < 25 kg/m2 or < 85th BMI percentile) Korean adolescents aged 13-16 years and subjected to 16S rRNA gene sequencing. Analysis of bacterial composition according to taxonomic rank (genus, family, and phylum) revealed marked differences in the Bacteroides and Prevotella populations in normal and obese samples (p < 0.005) at the genus and family levels; however, there was no difference in the Firmicutes-to-Bacteroidetes (F/B) ratio between normal and obese adolescents samples at the phylum level (F/B normal = 0.50 ± 0.53; F/B obese = 0.56 ± 0.86; p = 0.384). Statistical analysis revealed a significant association between the compositions of several bacterial taxa and child obesity. Among these, Bacteroides and Prevotella showed the most significant association with BMI (p < 0.0001 and 0.0001, respectively). We also found that the composition of Bacteroides was negatively associated with triglycerides (TG), total cholesterol, and high-sensitive C-reactive protein (hs-crp) (p = 0.0049, 0.0023, and 0.0038, respectively) levels, whereas that of Prevotella was positively associated with TG and hs-crp levels (p = 0.0394 and 0.0150, respectively). We then applied the association rule mining algorithm to generate "rules" to identify the association between the populations of multiple bacterial taxa and obesity; these rules were able to discriminate obese from normal states. Therefore, the present study describes a systemic approach to identify the association between bacterial populations in the gut and childhood obesity.

The feasibility study of non-invasive fetal trisomy 18 and 21 detection with semiconductor sequencing platform.

OBJECTIVE: Recent non-invasive prenatal testing (NIPT) technologies are based on next-generation sequencing (NGS). NGS allows rapid and effective clinical diagnoses to be determined with two common sequencing systems: Illumina and Ion Torrent platforms. The majority of NIPT technology is associated with Illumina platform. We investigated whether fetal trisomy 18 and 21 were sensitively and specifically detectable by semiconductor sequencer: Ion Proton. METHODS: From March 2012 to October 2013, we enrolled 155 pregnant women with fetuses who were diagnosed as high risk of fetal defects at Xiamen Maternal & Child Health Care Hospital (Xiamen, Fujian, China). Adapter-ligated DNA libraries were analyzed by the Ion Proton™ System (Life Technologies, Grand Island, NY, USA) with an average 0.3× sequencing coverage per nucleotide. Average total raw reads per sample was 6.5 million and mean rate of uniquely mapped reads was 59.0%. The results of this study were derived from BWA mapping. Z-score was used for fetal trisomy 18 and 21 detection. RESULTS: Interactive dot diagrams showed the minimal z-score values to discriminate negative versus positive cases of fetal trisomy 18 and 21. For fetal trisomy 18, the minimal z-score value of 2.459 showed 100% positive predictive and negative predictive values. The minimal z-score of 2.566 was used to classify negative versus positive cases of fetal trisomy 21. CONCLUSION: These results provide the evidence that fetal trisomy 18 and 21 detection can be performed with semiconductor sequencer. Our data also suggest that a prospective study should be performed with a larger cohort of clinically diverse obstetrics patients.